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The fate of our universe is an unceasing topic of cosmology and the human being. The discovery 
of the current accelerated expansion of the universe significantly changed our view of the fate of the 
universe. Recently, some interesting scenarios concerning the fate of the universe attracted much 
£U attention in the community, namely the so-called "Little Rip" and "Pseudo-Rip" . It is worth noting 

that all the Big Rip, Little Rip and Pseudo-Rip arise from the assumption that the dark energy 
density p(a) is monotonically increasing. In the present work, we are interested to investigate what 
will happen if this assumption is broken, and then propose a so-called "Quasi-Rip" scenario, which 
£\j " is driven by a type of quintom dark energy. In this work, we consider an explicit model of Quasi-Rip 

in details. We show that Quasi-Rip has an unique feature different from Big Rip, Little Rip and 
Pseudo-Rip. Our universe has a chance to be rebuilt in the ash after the terrible rip. This might 
^ | be the last hope in the "hopeless" rip. 
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I. INTRODUCTION 

Since its discovery in 1998, the current accelerated expansion ol our universe [l[ has been one ol 
the most active fields in modern cosmology. As is well known, it could be due to an unknown energy 
component (dark energy) or a modification to general relativity (modified gravity) [H, 0. Before 1998, 
it is commonly believed that in the future our universe will either expand forever or contract again 
into a final Big Crunch. However, the discovery of the current accelerated expansion of the universe 
significantly changed our view of the fate of the universe. In fact, many novel possibilities are under the 
active consideration in the community nowadays. 

Today, there are many dark energy candidates in the market. Among them, the famous phantom dark 
energy [3, 0| is very interesting. Its equation-of-state parameter (EoS) is smaller than —1. Although 
phantom dark energy is consistent with the current observational data [3, Q , it violates all the energy 
conditions. One of the consequences is that our universe will encounter a singularity at a finite time, 
namely the so-called Big Rip 4]. At this singularity, the scale factor a, energy density and pressure of 
our universe are all divergent. In fact, besides the traditional Big Bang, Big Crunch, and the Big Rip, 
many novel singularities have been considered in the literature, such as Sudden singularities, Generalized 
sudden singularities, Quiescent singularities, Big Boost, Big Brake, Big Freeze, w singularities, Inacces- 
sible singularities, Directional singularities (see e.g. 0, [a] and references therein for some brief reviews). 
These singularities arise at the price of violation of one or several energy conditions. As is well known, 
in Q (see also e.g. [1,0,0]) the future singularities have been classified into four types, namely 

• Type I (Big Rip): a — > oo, p — > oo, H — > oo, \p\ — > oo, when t — > t s < oc ; 

• Type II (Sudden singularity): a — > a s , p — > p s , H — > H s , \p\ — > oo, H — > oo, when t — > t s < oo ; 

• Type III (Big Freeze): a — ► a s , p — > oo, H — > oo, \p\ — ► oo, when t — > t s < oo ; 

• Type IV (Generalized sudden singularity): a — > a s , p — > p s , H H s , \p\ p s , H — > H s , and 
higher derivatives of H diverge, when t — > t s < oo , 

where t s , a s , p s , H s , p s , H s are all finite constants (a s ^ 0); p and p are energy density and pressure 
respectively; H = a/ a is the Hubble parameter, and a dot denotes a derivative with respect to the cosmic 
time t. In fact, the above four types can include almost all known future singularities. 

Of course, singularities usually are not desirable in physics. Therefore, other possible fates of our 
universe are also considered in the literature, such as the cyclic/oscillatory cosmology. Recently, some 
interesting scenarios concerning the fate of the universe attracted much attention in the community, 
namely the so-called "Little Rip" Q and "Pseudo-Rip" [l(|. In (9MTT| . their authors showed that if the 
cosmic energy density will remain constant or monotonically increase in the future, then all the possible 
fates of our universe can be divided into four categories based on the time asymptotics of the Hubble 
parameter H(t) [io| . namely 



Big Rip: H (<) — > oo, when t — > t r i p < oo ; 
Little Rip: H (t) — > oo, when t oo ; 
Cosmological Constant: H (t) — const. ; 
Pseudo-Rip: H(t) — > Hoo < oo, when t — > oo , 



where is a constant. Obviously, the Big Rip singularity is not the only fate of our universe with the 
phantom-like dark energy. Both Little Rip and Pseudo-Rip are non-singular, and hence fall outside of 
the four categories in Q. Similar to the Big Rip, the Little Rip dissociates all bound structures, but the 
strength of dark energy is not enough to rip apart spacetime (unlike the Big Rip) 0, [l(| • On the other 
hand, the Pseudo-Rip dissociates the bound structures which are held together by a binding force at or 
below a particular threshold, and hence it is possible that only some bound structures are dissociated 
while the others are not dissociated (depending on the model parameters) 10|. In fact, the Little Rip is 
an intermediate case between the cosmological constant and the Big Rip [9 1 , while the Pseudo-Rip is an 
intermediate case between the cosmological constant and the Little Rip jlO | . 
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It is worth noting that all the Big Rip, Little Rip and Pseudo-Rip arise from the assumption that 
the dark energy density p(a) is monotonically increasing 

Gum 

, i.e., the dark energy is phantom-like (its 
EoS w < —1). In the present work, we are interested to investigate what will happen if this assumption 
is broken. Obviously, in the case of the dark energy density p{a) is monotonically decreasing (i.e., the 
dark energy is quintessence-like with a EoS w > —1), no rip will happen and no bound structures will 
be dissociated. On the other hand, in the case of the dark energy density p(a) monotonically decreases 
(namely w > — 1) in the first stage and then monotonically increases (namely w < — 1) in the second stage 
(this is the case of the so-called "quintom A" dark energy in terminology of e.g. [l2| and references therein), 
the fate of our universe is the Big Rip, which is trivial in some sense. The third case is that the dark 
energy density p{a) monotonically increases (namely w < —1) in the first stage and then monotonically 
decreases (namely w > —1) in the second stage (this is the case of the so-called "quintom B" dark energy 
in terminology of e.g. [l2| and references therein). It can be expected that in the first stage some or 
all bound structures will be dissociated (similar to the case of Pseudo-Rip) , but then the disintegration 
process will stop, and the already disintegrated structures have the possibility to be recombined in the 
second stage. We dub it "Quasi- Rip" , which is the subject of the present work (here we temporarily do 
not consider the case of oscillatory quintom dark energy). Since the quintom-like dark energy [13j (whose 
EoS can cross the so-called phantom divide w = —1) is slightly favored by the observational data [lj (see 
e -g- UM f° r a comprehensive review), we note that the Quasi-Rip is well-motivated in fact. 

This paper is organized as followings. In SecHH we discuss the disintegration of bound structures. In 
Sec. IIIII we present an explicit model of Quasi-Rip. We constrain this model with the observational data, 
and then clearly show the Quasi-Rip in this model. In Sec. IIV[ some concluding remarks are given. 



II. THE DISINTEGRATION OF BOUND STRUCTURES 



It is useful to introduce the concept of "inertial force" when we discuss the disintegration of bound 
structures. According to [13, HH, H3 > m a flat Friedmann- Robertson- Walker (FRW) universe dominated 
by dark energy, (motivated by the well-known Newton's second law) the inertial force Fi„ ert on a mass 
m as seen by a gravitational source separated by a comoving distance ro is given by 

_ a ( ■ tt2 \ AttG , A-kG ( dp\ . 

F inert = mr o- = mr + H J = -mr — — (p + 3p) = mr — — I 2p + a— 1 , (1) 

in which we have used the energy conservation equation p + 3H(p+p) = 0. A bound structure dissociates 
when the inertial force Fi nert (dominated by dark energy) is equal to the force Fb ounc i holding together 
this bound structure [I(J El- F° r the gravitationally bound structure whose mass is M, it dissociates 
when 



a Mm 2 

Fmert = mr - = F bound = G — — = mr U} Q , (2) 
a r n 



namely 



or equivalently, 



'd 2 GM 
- =UJ o = — — 



a 



(3) 



'o 



1 ( n dp\ 1 / dp \ 2H^ 2 GM 

f(a) = - [2p + a-f)=- [2p+-^- )=— °- (4) 

po \ da J po \ din a J il rg 

where wo is the angular velocity; po is the present density of dark energy; f2 = (8nGpo)/ (3i?o) is the 
present fractional density of dark energy. Note that in this paper, the subscript "0" usually indicates 
the present value of corresponding quantity, and we have set do = 1. In fact, one can check that the 
results Eqs. ([3]) and (Q| are coincident with the ones in the case of phantom dark energy whose EoS w is 
a constant 0] (we refer to e.g. [TJ] for its detailed derivation; note that the result of Q 1& invalid for a 
non-constant w, while our results Eqs. ^ and Q are still valid). 
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However, if the bound structure is massive enough to significantly affect the local spacetime metric, it 
is not accurate to express Fi nert in terms of FRW metric 0, [13] ■ A more accurate method was presented 
in [TJ]. Following in the Newtonian limit, the interpolating metric of the local spacetime is given by 



ds 2 = 



! 2GM 

a(t)i] 



dt 2 - a 2 (t) [drj 2 + rf (d6 2 + sin 2 6dip 2 )] , (5) 



where rj is the comoving radial coordinate. The radial equation of motion for a test particle in the 
Newtonian limit reads 

r = -r+ — — , (6 

a r A 

where L = r 2 (p — r 2 ui = const. = TqWo is the constant angular momentum per unit mass. Noting that 
uj 2 = GM/ Tq , Eq. © can be recast as 



r = — r 



a GMr GM dV c 



cir 



(7) 



a r 3 r 2 dr 

where the time-dependent effective potential is given by 

18, GMrp GM 
VcS = -2a r+ — 2 (8) 

The bound structure dissociates when the minimum of the time-dependent effective potential (including 
the centrifugal term) disappears. Note that the corresponding radius r v at which the effective potential 
reaches its minimum is determined by 

^r* mjn - GMr Vmm + GMr Q = . (9) 

One can find (with the help of e.g. Mathematica) that Eq. (O has a real solution only for 

a . 27 GM 
r o 

Therefore, the minimum of the time-dependent effective potential disappears when 

d 27 GM (n) 



a 256 7-q 



which is also the time when the bound structure dissociates. Similarly, we can also recast Eq. (|TT|) as 



1 ( n dp\ 1 / dp \ 27 H7 2 GM 

{a ) = — (2p + a-f )=- [2p+-J- =— -A 3-. 12 

Po \ da J po \ d In a J 128 Hq 



As is shown in [14| , the accurate disintegration time obtained here is usually earlier than the qualitative 
one obtained in [4[ (this point can be seen also by comparing Eqs. (fTTj) . (fT2j) with Eqs. ©, (j4])). It is 
worth noting that most of the results in [bij are valid for dark energy whose EoS w is a constant, while 
our results obtained here are also valid for dynamical dark energy whose EoS w is not a constant. 

In the present work, we consider five bound structures, namely the Coma Cluster (M = 6 x 10 48 g, 
r = 9 x 10 24 cm), the Milky Way galaxy (M = 2 x 10 45 g, r = 5 x 10 22 cm), the Solar System 
(M = 2 x 10 33 g, r = 7 x 10 15 cm), the Earth (M = 6 x 10 27 g, r = 6.4 x 10 8 cm), and the Hydrogen 
atom. Note that the first three bound structures are suitable for Eqs. (ITT]) and (jT^J) , while the Earth is 
suitable for Eqs. ^ and (jU), since its surface material has no centrifugal term in the effective potential. 
Finally, the Hydrogen atom is also suitable for Eqs. ([3]) and (j4|), but in which the term GM should be 
replaced by a new term g 2 /(47reow e ), since Fb oun d is the electromagnetic force in this case. Note that 
q e = 1.6 x 10" 19 C, m e = 0.511 MeV (lMeV = 1.7827 x lCT 27 g), 1/(4tte ) = 9 x 10 18 gem 3 C" 2 sec" 2 , 
and r = 5.3 x 10~ 9 cm for the Hydrogen atom. In Eqs. (UJ) and (fT2"j). G — 6.672 x 10~ 6 cm 3 g _1 sec -2 , 
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H^ 1 = 3.0856 x 10 17 /i _1 sec (here h is the Hubble constant H in units of 100 km/sec/Mpc). Note that 
£7o and h will be determined by the observational data (see below). 

Since the current observational data are in the epoch a < 1, we cannot ignore the contribution from 
pressureless matter in this stage, although it can be safely ignored in the above discussions when the 
universe is dominated by dark energy. From the Friedmann equation 3H 2 = 8irGptot = StvG(p + p m ) and 
pm = Pmoa~ 3 , we have 



e 2 = = n — + n m0 a- 3 = n a £. + (i- n ) a ~ 3 , (13) 

Hq Po Po 



where £l m o = (^Gp m a) / (3Hq) is the present fractional density of pressureless matter. If p(a) is given, 
we can use Eq. (|13p to constrain the model parameters with the current observational data. On the 
other hand, from Eqs. (|12[) or one can find the corresponding scale factor a* at which the bound 
structure dissociates. Then, noting H = a/a, we can evaluate the disintegration time measuring from 
today (a — 1), namely 

da . f lna ' dlna 



where H^ 1 = 9.7776 h' 1 Gyr. 



Hn 1 / = H^ 1 / —, r , (14) 



III. AN EXPLICIT MODEL OF QUASI-RIP 



Here, we consider an explicit model of Quasi-Rip. In some sense, constructing a suitable model is a 
smart task. There are various ways to this end. For examples, one can specify the scale factor a as a 
function of the cosmic time t (see e.g. [EES 13); the pressure p as a function of the energy density p (see 
e.g. 0, H3, EH), the energy density p as a function of the scale factor a (see e.g. [j, fTol|V or the Hubble 
parameter if as a function of the cosmic time t (see e.g. [HI). It is shown in [9Hll| that these ways are 
equivalent in fact. In our case, for convenience, we choose to specify the e nerg y density p as a function 
of the scale factor a, similar to the cases of Little Rip [9j and Pseudo-Rip [10] . As mentioned in Sec. HI 
here our task is to construct an explicit function p{a), which monotonically increases (namely w < — 1) 
in the first stage and then monotonically decreases (namely w > —1) in the second stage. A naive idea 
is to use a piecewise function, which is phantom-like in the first sector (namely p(a) oc a~ 3 ( 1+Wl ' with a 
constant EoS w\ < —1) and is quintessence-like in the second sector (namely p(a) oc a^ 3t - 1+W2 ' ) with a 
constant EoS W2 > —1). Then, let us refine this naive idea with a smooth function p(a) cx o^', in which 
p(a < at) > and p(a > a t ) < 0, where at is the transition point. Noting that it is more convenient to 
use In a as the variable in Eqs. (fT12j). (U), and ([14")) , as well as the fact a — e lna , it is appropriate to write 
p as a function of In a instead. Therefore, the simplest function p(a) can be given by 

p(a) - Pa a a -t }lna = Po e lna i a ~P^ a ) , (15) 

where a and /3 are both constants. Let us have a closer observation. By requiring dp/dhia — 0, we find 
that its extremum locates at lna = a/(2/3). To ensure this extremum is a maximum, d 2 p/d(lna) 2 should 
be negative here, therefore /3 > is required. 
Substituting Eq. (fTS")) into Eq. (fTB"]). we have 

nQa «-pina + q_ q o j a -3 = Q Qe lna (a-/31na) + ^ _ Q ^ £ -3 In a 

no(i + z)-«~^+*) + (i _ n ) (i + z f , (i6) 

where z — 1/a — lis the redshift. There are three free model parameters, namely flo, a and /3. With this 
E(z), following the methodology in e.g. [HrHlj], we can constrain this model with the latest Union2.1 
Type la supernovae (SNIa) dataset 22] which consists of 580 SNIa. In fact, we find that the best fit 
has Xmin = 562.224, and the best-fit parameters are fl Q = 0.718859, a = 0.0294816 and (3 = 0.0002525, 
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whereas the corresponding h = 0.698862. In Fig. [TJ we present the corresponding 68.3% and 95.4% 
confidence level (C.L.) contours in the a — (3 parameter space. Although their best-fit values are small, 
a and (3 can still be consistent with the observational data in the large 2a C.L. region, i.e., — 1.2 < a < 0.5 
and < /3< 3.2. 

3.5 

3 
2.5 

2 
1.5 

1 

0.5 



-1 -0.5 0.5 

a 

FIG. 1: Fhe 68.3% and 95.4% confidence level contours in the a — /3 parameter space. The best-fit parameters 
are also indicated by a black solid point. 

Next, we show the Quasi- Rip in this model clearly. Substituting Eq. (|15p into the reduced inertial force 
/(a) defined in Eqs. (IT2"|) and (U]), we have 

/(a) = -^-(2 + a-2^1na) = (2 + a - 2/3 In a) a a ' plna = (2 + a - 2/31na) e lna (a ^ lna) . (17) 
Po 

In Fig. [21 we plot ln(p/po) and In / as functions of In a for various model parameters. For convenience, 
we fix j3 = 0.0001 which is at the same order of its best-fit value, while the various values of a are taken 
from the la C.L. region in Fig. Q] From the left panel of Fig.O it is easy to see that our p(a) in Eq. (fT5")) 
is the desirable one, which monotonically increases in the first stage and then monotonically decreases 
in the second stage. It is not surprising that the plot of In / is very similar to the one of ln(p//?o), since 
their difference is given by In (2 + a — 2/3 In a) which is a small quantity in fact (cf. Eq. (fTT)) ). Depending 
on the model parameters a and /3, the inertial force can successively disintegrate the Coma Cluster, the 
Milky Way galaxy, the Solar System, the Earth and the Hydrogen atom. In different cases, some bound 
structures will be disintegrated while the other bound structures will not. If none of the bound structures 
is disintegrated, the corresponding Quasi-Rip is a failed rip. From the right panel of Fig. [21 we find that 
the innermost one is a failed rip, the outermost one is the strongest Quasi-Rip which can disintegrate all 
the five bound structures, and the others are moderate ones which can only disintegrate one or several 
of the five bound structures. From Eqs. (1121) or (|4]), we can find the corresponding In a* at which the 
bound structure dissociates. Then, we can evaluate the disintegration time measuring from today i* — tg 
by using Eq. ([Tfj) with E'(lna) given in Eq. ([TBI) , and present the results in Tab.[TJ For a fixed /3, a bound 
structure dissociates earlier with a larger a, since the corresponding inertial force is stronger. In the case 
of a — 0.27, the difference between the disintegration time of Earth and Solar System is 4349 years, while 
the difference between the disintegration time of Hydrogen atom and Earth is only 15 years. Finally, it is 
easy to see that the most distinct feature of Quasi-Rip is that the inertial force monotonically decreases 
in the second stage. Eventually, it will become lower than all the thresholds to disintegrate the bound 
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structure. Therefore, the already disintegrated structures have the possibility to be recombined in the 
second stage. This is the unique feature of Quasi-Rip different from Big Rip, Little Rip and Pseudo-Rip. 
Our universe has a chance to be rebuilt in the ash after the terrible rip. This might be the last hope in 
the "hopeless" rip. 




FIG. 2: ln(p/po) and In/ as functions of lna for various model parameters. The curves from the innermost to 
the outermost are for a = 0.04, a = 0.064, a — 0.1, a = 0.15, a = 0.2 and a = 0.27, respectively, while /3 is fixed 
to be 0.0001. The dashed horizontal lines in the right panel indicate the corresponding thresholds to disintegrate 
the Coma Cluster, the Milky Way galaxy, the Solar System, the Earth and the Hydrogen atom (from bottom to 
top). Note that fio and h are taken to be their best-fit values. 



Bound structure 


q = 0.04 a = 0.064 a = 0.1 


a = 0.15 


a = 0.2 


a = 0.27 


Coma Cluster 


563.154 333.555 


216.991 


161.28 


118.715 


Milky Way 


344.035 


223.126 


165.75 


122.005 


Solar System 




223.281 


165.857985136 


122.083461362 


Earth 






165.857991583 


122.083465711 


Hydrogen atom 








122.083465696 



TABLE I: The disintegration time measuring from today t* — to (in units of Gyr) for models with various a, while 
j3 is fixed to be 0.0001, and Slo and h are taken to be their best-fit values. 



IV. CONCLUDING REMARKS 



The fate of our universe is an unceasing topic of cosmology and the human being. The discovery of the 
current accelerated expansion of the universe significantly changed our view of the fate of the universe. 
Recently, some interesting scenarios concerning the fate of the universe attracted much attention in the 
community, namely the so-called "Little Rip" and "Pseudo-Rip" . It is worth noting that all the Big Rip, 
Little Rip and Pseudo-Rip arise from the assumption that the dark energy density p{a) is monotonically 
increasing. In the present work, we are interested to investigate what will happen if this assumption is 
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broken, and then propose a so-called "Quasi-Rip" scenario, which is driven by a type of quintom dark 
energy. In this work, we consider an explicit model of Quasi-Rip in details. We show that Quasi-Rip has 
an unique feature different from Big Rip, Little Rip and Pseudo-Rip. Our universe has a chance to be 
rebuilt in the ash after the terrible rip. This might be the last hope in the "hopeless" rip. 

Some remarks are in order. Firstly, as is shown in Sec. HID our Quasi-Rip model is well consistent 
with the current observational data. However, even in the la C.L. region of a — (3 parameter space, the 
future behavior of our universe can be different enough (depending on the particular model parameters a 
and 0). In fact, as is well known, the current observational data can be consistent with all the phantom- 
like, quintessence-like and quintom-like dark energy models. Therefore, the current observational data 
cannot tightly tell what is the true fate of our universe. Most of the possibilities (including Big Rip, Little 
Rip, Pseudo-Rip, Quasi-Rip, dc Sitter expansion, other future singularities and so on) are still living. 

Secondly, the explicit model of Quasi-Rip considered in the present work is the simplest case. One can 
construct other more complicate p(a) to implement the Quasi-Rip. For example, one might construct 
an EoS w(a) as a function of scale factor a, which is smaller than —1 when a < at and is larger than 
— 1 when a > a t . Then the corresponding p(a) can be found from the energy conservation equation 
p + 3Hp(l + w(a)) = 0. Of course, other smart methods to construct the desirable p(a) are waiting us. 

Thirdly, as mentioned in in its second Pseudo-Rip model, the reduced inertial force f(a) can also 
have a peak, similar to our Quasi-Rip model. However, we note that in the Pseudo-Rip model, after the 
peak, the reduced inertial force /(a) — » const, which is still higher than the corresponding threshold to 
disintegrate the bound structure. Therefore, the already disintegrated structures have no possibility to 
be recombined in the Pseudo-Rip models. On the contrary, in the Quasi-Rip models, the reduced inertial 
force f(a) monotonically decreases in the second stage. Eventually, it will become lower than all the 
thresholds to disintegrate the bound structure. Therefore, the already disintegrated structures have the 
possibility to be recombined in the second stage. 

Fourthly, it is well known that phantom is unstable at quantum level and hence the perturbations 
grow large. Noting that in the present work our discussions are at classical level instead, this problem 
could be set aside. In fact, the quantum stability of a phantom phase has been considered in [24] ■ The 
authors of 24] studied the perturbations in the quantum-corrected effective field equation at one and 
two-loop order, and they found that the system is stable. On the other hand, it is claimed in (2f| that 
scalar perturbations can grow during a phantom phase if EoS w < —5/3. However, from Eq. (|15[) and 
Fig. [TJ it is easy to see that the corresponding w is only slightly smaller than —1 in the phantom phase 
for our particular a and (3 (n.b. Fig. [T|). So, w > —5/3 instead and hence our Quasi-Rip model can avoid 
the problem raised in [251 ]. Further, in fact the dark energy considered in this work is not necessary to 
be a scalar field. It could even be an effective dark energy from modified gravity, namely the so-called 
"geometric dark energy" . So, it might avoid the corresponding problems in the phantom phase. 

Finally, in the present work, we consider only the quintom dark energy which crosses the phantom 
divide w = — 1 once. In fact, we can also consider the quintom dark energy which can cross the phantom 
divide for many times. The most attractive Quasi-Rip model might be the one driven by the oscillatory 
quintom dark energy. In this oscillatory Quasi-Rip model, our universe will be destroyed and then be 
rebuilt again and again. 
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